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Abstract It is easy to sum chain-free self-energy rainbows, to obtain contributions to 
anomalous dimensions. It is also easy to resum rainbow-free self-energy chains. Taming 
the combinatoric explosion of all possible nestings and chainings of a primitive self-energy 
divergence is a much more demanding problem. We solve it in terms of the coproduct 
A, antipode S, and grading operator Y of the Hopf algebra of undecorated rooted trees. 
The vital operator is S*Y, with a star product effected by A. We perform 30-loop Pade- 
Borel resummation of 463 020 146 037 416 130 934 BPHZ subtractions in Yukawa theory, at 
spacetime dimension d = 4, and in a trivalent scalar theory, at d — 6, encountering residues 
of S-kY that involve primes with up to 60 digits. Even with a very large Yukawa coupling, 
g = 30, the precision of resummation is remarkable; a 31-loop calculation suggests that it 
is of order 10~ 8 . 
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1 Introduction 



In this work we develop the Hopf algebra of renormalization |T|, |2|, |||, |4], || to progress 
beyond the rainbow |], 0] and chain [§, ||] approximations for anomalous dimensions. 

Summing rainbows: In d dimensions, the massless scalar one-loop integral with 
propagators to the powers a, (3 is 

G(a, (3; d) := g(a)g((3)g(d - a - (3); g(a) := T(d/2 - a)/T(a) (1) 

Now consider the interaction g<jya<f>, with a neutral scalar particle a coupled to a charged 
scalar 0, in the critical dimension, d c = 6. To find the anomalous field dimension 7 of </>, 
in the rainbow approximation of 0, one solves the consistency condition 

1 = aG(l, 1 + 7; 6) = — a (2) 

717- l)(7-2)(7-3) 

which ensures that the coupling a := g 2 / (4:Ti) dc ^ 2 cancels the insertion of the anomalous 
self energy. The perturbative solution of the resulting quartic is easily found: 



5 + 4a/1 + a a _ a 2 
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7rainbow = " ~ = - - + 1 1 — - 206— H (3) 



6 6 3 6 5 

Resumming chains: At the other extreme, one may easily perform the Borel re- 
summation of chains of self-energy insertions, within a single rainbow. Suppose that the 
self energy p 2 ~E(a,p 2 / /j, 2 ) is renormalized in the momentum scheme, and hence vanishes 
at p 2 = fi 2 . The renormalized massless propagator is D = l/(p 2 — p 2 £). Then (§) is the 
rainbow approximation for <9£/<91og(/i 2 ) at p 2 = jj 2 . Following the methods of 0, one 
finds that the corresponding asymptotic series for chains is Borel resummable: 

Jo (x + l)(x + 2)(x + 3) 6 6 6 6 b 

which differs from the rainbow approximation at 3 loops, with 206 in @ coming from the 
triple rainbow, while 170 in (Q) comes from a chain of two self energies inside a third. 

Hopf algebra: We shall progress beyond the rainbow and chain approximations by 
including all possible nestings and chainings of the one-loop self-energy divergence. In 
other words, we consider the full Hopf algebra of undecorated rooted trees, established 
in [|l| and implemented in @]. Two figures suffice to exhibit the class of diagrams con- 
sidered, and their divergence structure. The first exhibits a 12-loop example, the second 
exhibits its divergence structure. Due to the fact that we combine chains and rainbows, 
we have a full tree structure the depth of the tree is larger than one, and there can 
be more than one edge attached to a vertex. 

There are 4 notable features of this analysis. 

1. We use the coproduct A to combine the antipode S and grading operator Y in a 
star product S *Y whose residue delivers the contribution of each rooted tree. 
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Figure 1: A 12-loop diagram based on a one-loop skeleton. 




2. We show that the rationality of rainbows [10 extends to the contribution of every 



undecorated rooted tree, as had been inferred from examples in [11 



3. We confirm that a recent analysis JT2 ] of dimensional regularization applies at both 
6, detecting poles of T functions that occur in even dimensions. 



d c = 4 and d c 



4. We obtain, to 30 loops, highly non-trivial alternating asymptotic series, which we 
resum, to high precision, by combining Pade [O and Borel H, M methods. 



2 Hopf-algebra method 

Let t be an undecorated rooted tree, denoting the divergence structure of a Feynman 
diagram. Then its coproduct is defined, recursively, by 

A(t)=t®e + id®B + (A(B-(t))) (5) 

where e is the empty tree, evaluating to unity, id is the identity map, 5_ removes the 
root, giving a product of trees in general, and B + is the inverse of -B_, combining products 
by restoring a common root. The recursion terminates with A(e) = e ® e and develops a 
highly non-trivial structure by the operation of the coproduct on products of trees 

a en***) = n fc A(t fc ) (6) 
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between each removal and restoration of a root. In Sweedler notation, it takes the form 

A(t) = Efc4 1} ® 4 2) = * ® e + e ® i + EfcO^ ® a? (7) 
with single trees on the right and, in general, products on the left. The prime in the 
second summation indicates the absence of the empty tree. The field-theoretic role of the 
coproduct is clear: on the left products of sub divergences are identified; on the right these 
shrink to points. Subtractions are effected by the antipode, defined by the recursion 

s(t) = -t- e^(4 1} )4 2) (8) 

for a non-empty tree, with S(Y\ k tk) = Yik S{t k ) for products and S(e) = e. 

Renormalization involves a twisted antipode, Sr. Let denote the Feynman map 
that assigns a dimensionally regularized bare value <f>(t) to the diagram whose divergence 
structure is labelled by the tree t. Then we apply the recursive definition Q 

S R (t) = -R (<j>{t) + E' fc ^(4 1) )0(4 2) )) (9) 

with a renormalization operator R that sets p 2 = [i 2 , in both the momentum and MS 
schemes, and in the MS scheme selects only the poles in e := (d c — d)/2. 

We can use the coproduct to combine operators. Suppose that 0\ and 2 operate on 
trees and their products. Then we define the star product 0\ * 2 by 

O l *0 2 {t) = Y jk 1 {a k l) )0 2 {a k 2) ) (10) 

with ordinary multiplication performed after 0\ operates on the left and 2 on the right 
of each term in the coproduct. By construction, S-kid annihilates everything except the 
empty tree, e. The presence of R makes Sr* <p finite and non-trivial. In particular, the 
renormalized Green function is simply 

r B (*) = limSW(*) (11) 

e^O 

whose evaluation was efficiently encoded in ||, using a few lines of computer algebra. 

Here we present a new - and vital - formula for efficiently computing the contribution 
of an undecorated tree to the anomalous dimension. It is simply 

7 (t)=lim^(S'*Y(t)) (12) 

where Y is the grading operator, with Y(t) = nt, for a tree with n nodes. In general, Y 
multiplies a product of trees by its total number of nodes. To see that this works, consider 
the terms in (0), in the momentum scheme, before taking the limit e — > 0. Each term has 
a momentum dependence (p 2 )™^ dc ^ 2 , where n is the number of loops (and hence nodes) 
of the tree on the right of the term in the Sweedler sum. If we multiply by ne, and then 
let e — > 0, we clearly obtain the derivative w.r.t. log(/i 2 /p 2 ). Setting p 2 = fi 2 we obtain the 
contribution to the anomalous dimension. Thus R plays no role and we may replace S#(t) 
by \im.R^ id Sr(J:) = <p(S(t)), where S is the canonical antipode. Multiplication by ne is 
achieved by £<f>(Y(t)) = ne<p(t) on the right of the coproduct, where Y acts only on single 
trees. Hence the abstract operator S * Y delivers the precise combination of products 
of trees whose bare evaluation as Feynman diagrams is guaranteed to have merely a 1/e 
singularity, with residue equal to the contribution to the anomalous dimension. Thus we 
entirely separate the combinatorics from the analysis. 
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3 Example 



By way of example, we show how the 3-loop expansions of result from fll2l) . The 



combinatorics are now clear. The analysis, at first sight, seems to entail the detailed 
properties of T functions. However, appearances can be misleading. 

In general, a dimensionally regularized bare value for a n-loop diagram, corresponding 
to the undecorated rooted tree t, is evaluated by the recursion |3J] 

0(t) = ^^n^) as) 

ne A fc A 

where hk are the branches originating from the root of t. It terminates with 0(e) = 1. For 
the scalar theory with d c = 6, the master function is 

r/ o a* , r G r(i - <j)r(i + <j)r(2 - e) 

Now the wonderful feature of ( |T2| ) is that it depends only on the derivatives of L(e, 5) 
w.r.t. 5 at e = 0. This reflects the fact that the anomalous dimension, unlike the Green 
function, is insensitive to the details of the regularization method. Thus we may, with 
huge savings in computation time, replace the master function by 

L( °' *> = (*-!)(*! 2)0- 3) = g**" = "5 + " < + °<*> < 15 > 

which establishes that the contribution of each rooted tree is rational. The residue of the 
anomalous dimension operator S * Y feels only the rational residues of T functions; it is 
blind to the zeta-valued derivatives that contribute to the renormalized Green function. 

Now that the analysis has been drastically simplified, we return to the combinatorics. 
The double rainbow, t 2) has coproduct A(t 2 ) = t 2 ®e + e®£2 +t% <8>ti where t\ is the single 
rainbow, with A(t x ) = t\®e + e®t\. The antipodes are S(ti) = —t\ and S(t 2 ) = —t 2 + t\. 
The star products are S*Y{tx) = t\ and S , *y(t 2 ) = 2t 2 —t\- Hence the contributions to 
the anomalous dimensions are the residues of L(0, e)/e and (L(0, 2e) — L(0, e))L(0, e)/s 2 , 
namely go = —a/6 and gigo = lla 2 /6 3 . 

Following this simple example, the reader should find it easy to determine the anoma- 
lous dimension contributions of the two rooted trees at 3 loops. For £3, the triple rainbow 
graph, S * Y delivers 3^3 — 3t]t 2 + tf, with residue g 2 gl + gfgo = —(85 + ll 2 )a 3 /6 5 , in 
agreement with (|3]). For the other diagram, t' 3 , with a double chain in a single rainbow, 
it delivers 3t' 3 — 4iit 2 + tf with residue 1g 2 g\ = —2 x 85a 3 /6 5 , in agreement with (|). The 
Borel resummation (ffi) of chains corresponds to the result n\g n gQ for a chain of n self 
energies, inside a single rainbow. Writing the anomalous dimension contribution of the 
full Hopf algebra as the asymptotic series 

7hopf - E (16) 

n>0 U 



4 



we find that G 3 = 3 x 85 + ll 2 = 376. 

In this paper, we undertake Pade-Borel resummation of the full Hopf series fllB]), to 
30 loops. We also resum 

Thopf ^ £ <5 tg: (it) 

n>0 Z 

for the anomalous dimension of a fermion field with a Yukawa interaction gipcrip, at d c = 4, 
whose rainbow approximation 

7rainbow = 1 ~ VT+O (18) 

was obtained in 0. At the other extreme, the Borel-resummed chain approximation 

/-°° exp(-2x/a)tfa 

7chain = -2 / — (19) 

Jo x + 2 

is easily obtained from the Yukawa generating function, L(0, 5) — a/(S — 2). 



4 Results to 30 loops 

At 4 loops, there are 5 undecorated Wick contractions, corresponding to 4 rooted trees, 
one of which has weight 2. For the scalar theory, at d c = 6, the tally is 

G 4 = 4890 + 4711 + 3595 + 3595 + 3450 = 20241 = 3 2 x 13 x 173 (20) 

Already this becomes tedious to compute by hand. Fortunately, the recursions (||@) of 
the coproduct and antipode make it sublimely easy to automate the procedure ([12]). 



At n loops, the number of relevant Wick contractions is the Catalan number C n _i, 
where C n := ^ T ( 2 „"). At 30 loops, there are C 29 = 1002 242 216 651368 contrac- 
tions. Symmetries reduce these to rooted trees, with weights determined recursively 
by W(t) = w{t)Y[ k W{bk) where bk are the branches obtained by removing the root of 
t. The symmetry factor of the root is w(t) = Ilj n y where rij is the number 

of branches of type j. The generating formula for R n , the number of rooted trees with 
n nodes, is |L4| I] n>0 -Rn^ n = ^rin>o(l — which expresses the fact that removal 



of roots from all trees with n nodes produces all products of trees with a total of n — 1 
nodes. This gives i? 30 = 354 426 847 597. The number of terms produced by applying the 
BPHZ procedure to a single tree with n nodes is 2". 

From these enumerations, one finds - with some trepidation - that computation to 
30 loops entails S n <30 2 n -R n = 463 020 146 037416 130 934 subtractions, each requiring 30 
terms in its Laurent expansion, with coefficients involving integers of O(10 60 ). Brute 
force would require processing of O(10 24 ) bits of data, which is far beyond anything 
contemplated by current computer science. The remedy is clear: recursion of coproduct 
and antipode, to compute the residues of the anomalous dimension operator S * Y. 

Each new coproduct or antipode refers to others with fewer loops. By storing these 
we easily progressed to 13 loops, extending the sequence G n to 
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1, 11, 376, 20241, 1427156, 121639250, 12007003824, 1337583507153, 
165328009728652, 22404009743110566, 3299256277254713760, 
524366465815117346250, 89448728780073829991976 



For G n , in the Yukawa case, we obtained the 13-loop sequence 

1, 1, 4, 27, 248, 2830, 38232, 593859, 10401712, 202601898 
4342263000, 101551822350, 2573779506192 

At this point, recursion of individual trees hit a ceiling imposed by memory limitations. 

Beyond 13 loops, we stored only the unique combination of terms that is needed 
at higher loops, namely the momentum-scheme renormalized self energy. Allocating 
750 megabytes of main memory to Reduce 3.7 |16|], the time to reach 30 loops was 8 hours. 



Of these, more than 2 hours were spent on garbage collection, indicating the combina- 
toric complexity. Results for the scalar and Yukawa theories are in Tables 1 and 2. They 
are highly non-trivial. Factorization of G 27 = 2 6 x 5 x 103 x 184892457645048836717 x 
69943104850621681268329469624581 needed significant use of Richard Crandall's elliptic 
curve routine [17 , while G29/24O is a 60-digit integer that is most probably prime. 



5 Pade-Borel resummation 



We combine Pade-approximant [13| and Borel-transformation §, || methods. From 



we obtain the pure chain contribution G^™ = (2 n + (2 n — l)3 n+1 )n! with, for example, 
^<cham = (8 + 7 x 81) x 6 = 3450 appearing in ( p0|) as the smallest contribution of the 
5 Wick contractions at 4 loops, while the pure rainbow contribution, 4711, is next to 
largest. This is far removed from the situation at large n, where the pure rainbow term is 
factorially smaller than the pure chain term. At large n, we combine C n -\ ~ 4 n ~ 1 /v / n 3 7r 
Wick contractions, some of which are of order G n ham , while some are far smaller. It 
is thus difficult to anticipate the large-n behaviour of G n . We adopted an empirical 
approach, finding that S n := 12 l ~ n G n /T{n + 2) varies little for n G [14,30], as shown in 
the final column of Table 1. In the Yukawa case of Table 2, we found little variation in 
S n :=2 1 -™G n /r(n + l/2). 

In the scalar case, at d c = 6, Pade-Borel resummation may be achieved by the Ansatz 



a 

7hopf « 



/ P{ax/3)e- x x 2 dx (21) 
Jo 



12 Jo 

where P(y) = 1 + 0(y) is a [M\N] Pade approximant, with numerator 1 + J2m=i c m,y m 
and denominator 1 + X^=i d n y n , chosen so as to reproduce the first M + N + l terms in the 
asymptotic series (0). We expect P(y) to have singularities only in the left half-plane. In 
particular, a pole near y = — 1 is expected, corresponding to the approximate constancy 
of S n in Table 1. We fitted the first 29 values of G n with a [14\14] Pade approximant 
P{y), finding a pole at y ~ —0.994. The other 13 poles have 3ly < —1. Moreover there is 
no zero with ^ty > 0. The test-value G30 is reproduced to a precision of 5 x 10~ 16 . 
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In the Yukawa case, at d c = 4, we made the Ansatz 

7ho P f « -4= [°° QiaxWe-'x 1 ' 2 dx; Q(y) := (22) 
V7T Jo 1 + y 



suggested by Table 2. Here we put in by hand the suspected pole at y — — 1. The [14\14] 
approximant to P{y) = 1 + 0(y) then has all its 14 poles at dty < — 1 and no zero with 
9fa/ > 0. The test-value G 30 is reproduced to a reassuring precision of 4 x 1CT 17 . 

Table 3 compares resummation of the full Hopf results (|1(||17]) with those from the far 
more restrictive chain and rainbow subsets. To test the precision of resummations (|21~| , |22"|) , 
we used the star product (12) to perform the 2.6 x 10 21 BPHZ subtractions that yield the 



exact 31-loop coefficients 

G?3i = 2 6 x 3 3 x 5 x 139 x 2957 x 22279 x 69318820356301 x 9602299922477621 

x 144927172127490232568467 (23) 
G 31 = 2 5 x 3 4 x 5 x 71 x 109 x 13224049649 x 473202021103152647613521 (24) 

No change in the final digits of Table 3 results from using these. At the prodigious Yukawa 
coupling g = 30, corresponding to a = (30/47r) 2 « 5.7, a [15\15] Pade approximant 
gives 7hopf ~ —1.85202761, differing by less than 1 part in 10 8 from the [14\14] result 
7ho P f ~ —1.85202762. It appears that resummation of undecorated rooted trees is under 
very good control, notwithstanding the combinatoric explosion apparent in (|2"3],|2"4l) . 



6 Conclusions 

As stated in the introduction, we achieved 4 goals. First, we found the Hopf-algebra 
construct (12) that delivers undecorated contributions to anomalous dimensions. Then 
we found that these are rational, with the T functions of (|I4"D contributing only their 
residues, via (0). Next, we exemplified the analysis of dimensional regularization in fl2|l , 



at two different critical dimensions, d c = 6 and d c = 4. The residues of a common set ([I]) 
of T functions determine both results. Finally, we obtained highly non-trivial results, 
from all combinations of rainbows and chains, to 30 loops. A priori, we had no idea how 
these would compare with the easily determined pure chain contributions. Tables 1 and 2 
suggest that at large n the full Hopf-algebra results exceed pure chains by factors that 
scale like n 2 2 n and n 1//2 2 n , respectively. Pade approximation gave 15-digit agreement with 
exact 30- loop results. In Table 3, we compare the Borel resummations (|2l|j22]) of the full 



Hopf algebra with the vastly simpler rainbow approximations (|3|,|T8"D and the still rather 
trivial chain approximations (^ JT9|) . Even at the very large Yukawa coupling g = 30 we 
claim 8-digit precision. Apart from large- Nf approximations [ I3fl , we know of no other 



large-coupling analysis of anomalous dimension contributions, at spacetime dimensions 
d > 4, that progresses beyond pure rainbows [H, [7| or pure chains 0, M. 

In conclusion: Hopf algebra tames the combinatorics of renormalization, by disentan- 
gling the iterative subtraction of primitive sub divergences from the analytical challenge 
of evaluating dimensionally regularized bare values for Feynman diagrams. Progress with 
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the analytic challenge shall require the expansion of skeleton graphs in the regularization 
parameter D — 4. After that, the Hopf algebra of decorated rooted trees provides the 
tool to take care of the combinatorial challenge of renormalization in general. Generaliza- 
tions of the methods here to cases where decorations are different, but still analytically 



trivial, are conceivable. The results in [|nj are of this form. In the present case, where 
the combinatoric explosion is ferocious, while the analysis is routine, the automation of 
renormalization by Hopf algebra is a joy. How else might one resum 2.6 x 10 21 BPHZ 
subtractions at 31 loops and achieve 8-digit precision at very strong coupling? 

Acknowledgements: This work was undertaken during the workshop Number Theory 
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Cartier, Werner Nahm, Ivan Todorov and Jean-Bernard Zuber. Past work with Alain 
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Table 1: Scalar coefficients in flT|), with S n := 12 l - n G n /T(n + 2) 



n 


G n 


S n 


14 


16301356287284530869810308 


0.1165 


15 


3161258841758986060906197536 


0.1177 


16 


650090787950164885954804021185 


0.1186 


17 


141326399508139818539694443090940 


0.1194 


18 


32389192708918228594003667471390750 


0.1200 


19 


7805642594117634874205145727265669184 


0.1205 


20 


1973552096478862083584247237907087008846 


0.1209 


21 


522399387732959889862436502331522596697560 


0.1212 


22 


144486332652501966354908665093390779463113660 


0.1215 


23 


41681362292986022786933211385817840822702468640 


0.1217 


24 


12520661507532542738174037622803485508817145773050 


0.1218 


25 


3910338928202486568787314743084879349561179264255736 


0.1220 


26 


1267891158800355844456289086726128521948839015617187260 


0.1221 


27 


426237156086127437403654947366849019736474802601497417920 


0.1221 


28 


148382376919675149120919349602375065827367635238832722748020 


0.1222 


29 


53428133467243180546330391126922442419952183999220340144106320 


0.1222 


30 


19876558632009586773182109989526780486481329823560105761256963720 


0.1222 
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Table 2: Yukawa coefficients in (JT7|), with S n := 2 1 - n G n /T{n + 1/2) 



n 


C n 




14 


70282204726396 


0.3715 


15 


2057490936366320 


0.3750 


16 


64291032462761955 


0.3780 


17 


2136017303903513184 


0.3806 


18 


75197869250518812754 


0.3828 


19 


2796475872605709079512 


0.3848 


20 


109549714522464120960474 


0.3865 


21 


4509302910783496963256400 


0.3880 


22 


194584224274515194731540740 


0.3894 


23 


8784041120771057847338352720 


0.3906 


24 


414032133398397494698579333710 


0.3917 


25 


20340342746544244143487152873888 


0.3928 


26 


1039819967521866936447997028508900 


0.3937 


27 


55230362672853506023203822058592752 


0.3946 


28 


3043750896574866226650924152479935036 


0.3953 


29 


173814476864493583374050720641310171808 


0.3961 


30 


10272611586206353744425870217572111879288 


0.3968 



Table 3: Comparison of chain, rainbow and full Hopf contributions 



a 




Ochain 




Orainbow 




7hopf 




Tchain 




^rainbow 




7hopf 


0.5 





0727579 





0731322 





0742476 





2245593 





2247449 





2278233 


1.0 





1301409 





1322419 





1373080 





4126913 





4142136 





4281423 


1.5 





1773375 





1825988 





1937609 





5765641 





5811388 





6118625 


2.0 





2172313 





2268615 





2455916 





7226572 





7320508 





7837372 
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3394353 
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1 
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3961033 
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2904639 
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3452079 
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4494897 


1 
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4694751 


1 
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1 
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1 
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1 
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1 
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